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ABSTRACT 



Comet P/2010 A2 LINEAR is an object on an asteroidal orbit within the inner Main Belt, therefore a good candidate for membership 
with the Main Belt Comet family. It was observed with several telescopes (ESO New Technology Telescope, La Silla, Chile; Gemini 
^ ^ ' North, Mauna Kea, Hawai'i; University of Hawai'i 2.2 m, Mauna Kea, Hawai'i) from 14 Jan. until 19 Feb. 2010 in order to characterize 

O ^. and monitor it and its very unusual dust tail, which appears almost fully detached from the nucleus; the head of the tail includes two 

I . narrow arcs forming a cross. No evolution was observed during the span of the observations. Observations obtained during the 

O ' Earth orbital plane crossing allowed an examination of the out-of-plane 3D structure of the tail. The immediate surroundings of the 

^ I nucleus were found dust-free, which allowed an estimate of the nucleus radius of 80-90 m, assuming an albedo p - 0.11 and a 

C/) , phase correction with G = 0.15 (values typical for S-type asteroids). A model of the thermal evolution indicates that such a small 

. nucleus could not maintain any ice content for more than a few million years on its current orbit, ruling out ice sublimation dust 

ejection mechanism. Rotational spin-up and electrostatic dust levitations were also rejected, leaving an impact with a smaller body 
as the favoured hypothesis. This is further supported by the analysis of the tail structure. Finston-Probstein dynamical dust modelling 
indicates the tail was produced by a single burst of dust emission. More advanced models (described in detail in a companion paper), 
independently indicate that this burst populated a hollow cone with a half-opening angle a ~ 40° and with an ejection velocity 
Vmax ~ 0.2 m s~\ where the small dust grains fill the observed tail, while the arcs are foreshortened sections of the burst cone. 
The dust grains in the tail are measured to have radii between a = 1-20 mm, with a diff'erential size distribution proportional to 
^-3.44+0.08 contained in the tail is estimated to at least 8x10^ kg, which would form a sphere of 40 m radius (with a density 

p = 3 000 kg m~^ and an albedo p = 0.11 typical of S-type asteroids). Analysing these results in the framework of crater physics, 
■ we conclude that a gravity-controlled crater would have grown up to ~ 100 m radius, i.e. comparable to the size of the body. The 

non-disruption of the body suggest this was an oblique impact. 

Key words. Comets: P/2010 A2 (LINEAR), Asteroids: P/2010 A2 (LINEAR), Techniques: image processing, photometric 

K> . 1 ■ Introduction envel ope and oxides in a magma ocean [Oenda & Ikomal (l2007l 

}^ l2008l) . or (///) delivery by volatile-rich planetesimals (asteroid s 

^ Habitability within our solar system is determined by the dis- ^^^^^^^ ^^^^^^ ^^^^^^ ^^^^ ^-^^ (iMorbidelli et al.l2000h . 

- - tribution of water and volatiles, yet the origin of this distribu- ^he first two processes probably cont ributed to, but ma y be un- 

tion IS currently a fundamental unresolved planetary science is- ^^^^ ^^^^^^^ of Earth's water (iMottl et al.l2007h . Within 

sue ( |Kasting&Catling| | 2003D . There are three leading scenarios ^^^^^^^ ^^^^^^^ ^^^-^^^ ^^-^ ^^^^ ^^^^^^ (MBCs), 

for the origin of terrestrial planetary water, including: (/) nebu- ^ ^^^^^ discovered class of "comets having sta ble orbits com- 

lar gaAadsorption on micron-sized dust grains inside the snow confined to the main asteroid belt" (iHsieh & JewittI 

line ( lMuralidharanetal. 1 [20081 the distance from the Sun out- ||^^ ^^^^^^^ ^ ^^^^^^^^ particular significance to the history 

side of which water the temperature is low enough for water ice terrestrial water and other important volatiles. 
to condense), (//) chemical reactions between an early hydrogen 

In addition to possibly P/2010 A2 (LINEAR), there are 



* Based on observations collected at the Gemini North Observatory, currently five known MBCs: 133P/Elst-Pizarro, 176P/LINEAR, 

Mauna Kea, Hawai'i, USA, program GN-2009B-DD-10 at the 238P/Read, C/2008 Rl (Garradd) and P/201Q R2 (La Sagra). 

European Southern Observatory, La Silla, Chile, program 184.C- The mam belt comets are defined (|Hsieh & Jewitt 2006) by (/) a 

1 143(A), and at the University of Hawai'i 2.2-m telescope, Mauna Kea, semi-major axis that is less than Jupiter's, (//) Tisserand parame- 

Hawai'i, USA. ters significantly greater than 3 — meaning they are dynamically 
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decoupled fr om Jupiter, like ordinary asteroids (IVaghil IT973I: 
lKresaklll98Ql) — and (///) mass loss with a cometary appearance. 
As comets in near-circular orbits within the asteroid belt, these 
objects likely still harbor n ebular water frozen out from beyond 
the primordial "snow line" (lEncrena 32008l:[G^d & Linl2007h 
of the young solar syste m. Dynamical simula tions suggest they 
probably formed in- situ (lHaghighipoud2009b at a different tem- 
perature from either comets or the asteroidal reservoir that has 
been sampled through our meteorite collections. Thus, under- 
standing their chemistry can provide unique insight into the dis- 
tribution of volatiles in the early stages of planet formation. 

P/2010 A2 (LINEA R) was discovered by the LINEAR 
project on 7 Jan., 2010 (iBirtwhistle et al.|[2010ah and was de- 
scribed as "a head less comet with a straig ht tail, and no cen- 
tral condensation" (IBirtwhistle et al.ll2010bl) . On 1 1-12 Jan., ob- 
servers at the WYNN 3.5-m and at the 2.5-m Nordic Optical 
Telescope reported that the object had a asteroidal-like body 
~ 150-200m in diameter which was co nnected to the tail by 
an unresolved light bridge (lGreenll2010h . Based on the orbital 
semi-major axis and on the Tisserand parameter, Tj = 3.6, 
Jewitt ( Green 2010) concluded that P/2010 A2 was a new MBC. 
With the smallest perihelion distance of any of the MBCs (q = 
2.29 AU), these elements suggested a membership in the Flora 
family. Flora is a large asteroid family of ~500 members which 
can be broken up into many sub-families or clans and which are 
broadly compositi onally consistent wit h space- weathered S- type 
asteroidal spectra (iFlorczak et al.l 19981) . Jewitt (lGreenl2010h and 
Licandro (iHaver etal.ll2010h sussested that the location of the 
nucleus outside the coma might be the consequence of an im- 
pact. 

Is P/2010 A2 a genuine, sublimating comet (whose activity 
was possibly triggered by an impact), or is the tail the signa- 
ture of an impact (or another alternative process), with a disper- 
sion of the dust ejecta but no sublimation? If it can be demon- 
strated that the object is a comet, it would indicate that at least 
one object in the inner asteroid belt still contains volatile ma- 
terial. The water snow line in the proto- solar ne bula is esti- 
mated to have been around 2.5 AU from the Sun ('Jones et alj 
19901) or possibly even as close as 1.5 AU ([Lecar et al. 2006; 
Machida& Ab3l2010l: iMin et al.ll20Tl1) . Initially parent bodies 
beyond the snow line are believed to be mixtures of water ice and 
silicates, which are then heated due to the radioactiv e decay of 
^^Al within 1-2 m illion years after nebular collapse (iKrot et al.l 
l2006h . iGrimm & M cSween (1989) find that once the water is 
liquid, it is consumed by hy dratio n reactions, prefer entially in 
the interior (ICohen & Cokeijl2000l: Iwilson et al.lll999h . possibly 
leaving ice in the outer layers. In the inner asteroid belt, it is 
likely that most of that water has bee n converted in hydr ated ma- 
terial observed on S-type asteroids (iRivkin et al.ll2002h . MBCs 
may be the frozen components of the outer edges of asteroid 
parent bodies that have survived the age of the solar system. 

Four of the MBCs are located in the outer main belt, where 
asteroids with hydrated material are less common, suggesting 
water ice could still exist. The orbital elements of the MBCs 
are listed in Table [T] and displayed in Fig.[T] 133P/Elst-Pizarro, 
the first MBC discovered, is a member of the Themis fam- 
ily; 176P/LINEAR was found in a survey tar geting object s 
with orbits similar to that of 133P/Elst-Pizarro (lHsiehll2009h . 
238P/Read, however, was discovered serendipitously (i.e. not 
within a MBC-dedicated program) on a similar orbit. A fourth 
MBC, P/2008 Rl Garradd, was discovered serendipitously in the 
central region of the main belt. The fifth MBC, P/2010 R2 La 
Sagra, was also discovered serendipitou sly with a semi-majo r 
axis similar to that of 133P/Elst-Pizarro (iMarsden et al.ll201Qh . 



Fig. 1. (not available on astroPh) The orbital elements of 
P/2010 A2 (green square to the left), of the 5 known 
MBCs (circles, 133P and 23 8P are indistinguishable) and of 
the numbered main belt asteroids (small dots). When avail- 
able, the proper elements were used (source: the Asteroid 
Dynamics Site, http://hamilton.dm.unipi.it/astdys/). The inclina- 
tion of P/2010 R2 La Sagra is off' scale; its a is marked by a 
triangle on the top of the plot; / = 21.39 . 



More rec ently, asteroid ( 596) Scheila presented a co met-like ap- 
pearance (lLarsonll201()h . lBodewits et al.l (l201lh and I Jewitt et al.l 
(201il) conclude however that this dust cloud was likely caused 
by the impact of a small asteroid on Scheila. The survival of an 
icy body in the inner asteroid belt would therefore give exciting 
constraints on the evolution of water in the belt. It would also 
raise fundamental questions on how to shield water ice in that 
area of the solar system in a rather small body. 

iMoreno et al.l (l2010l) presented observations obtained with 
the GTC, WHT and NOT on La Palma, which they mod- 
ele d with an extended period of water-driven cometary activ- 
ity. iJewitt et aP (l2010l) acquired a series of HST images; from 
the orientation and geometry of the tail, they favor the disrup- 
tion of an ast eroid (either by collisio n or spin up) in Feb.-Mar. 
2009. Finally, [Snodgrass et al.l (l2010l) secured observation from 
the Rosetta spacecraft. Thanks to the position of the space probe, 
they could observe the object with a very diff'erent geometry (but 
with a much more modest resolution), from which they con- 
cluded that the observed dust tail was caused by an impact. 

In order to investigate the process that generated the ob- 
served dust tail, we acquired deep images of P/2010 A2 over 
various epochs. The observations are presented in Section[2j The 
analysis and modelling of the nucleus and dust are described in 
Section [3l In Section |4l we discuss th e conclusions and s umma- 
rize the results. A companion paper (iKleyna et al.ll201 ll) is de- 
voted to the details of the dust n iodels developed for th is study, 
and a second follow-up article (iHermalyn et al.|[20Tl1) focuses 
on the characteristics of the impact process and of the nucleus 
based on the interpretation of the dust as an impact plume. 



2. Observations 

The telescopes and instruments used are described below. The 
epoch, geometric circumstances and a log of the observations 
are listed in Table O 

In all cases, the observations were acquired as series of fairly 
short individual exposures obtained while tracking at the non- 
sidereal motion rates of the comet, and off'setting the telescope 
position between exposures. 

2.1. Telescopes and Instruments 
2.1.1. New Technology Telescope 

The observations were performed on the ESO 3.56m New 
Technology Telescope (NTT) on La Silla, with the ESO Faint 
Object Spectrograph and Camera (v . 2) ins trument (EF0SC2) 
(iBuzzoni et al ] |1984I: ISnodgrass et ani2008h . through Bessel B, 
y, R, and Gunn / filters, using the ESO#40 detector, a 2kx2k 
thinned, UV-ffooded Loral/Lesser CCD, which was read in a 2x2 
bin mode resulting in Off 24 pixels, and a 4f 1 field of view. 
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Fig. 2. P/2010 A2, images from UT 14.2 (a), 16.2 (b), 17.2 (c), 
18.2 (d) Jan. 2010 using the NTT, UT 19.5 Jan. 2010 using GN 
(e). The Hnear gray scale covers the range of (0-3) xlO"^A/, a 
proxy for the amount of dust present (see Section [3.3.3b . Each 
panel includes a 40x30'' zoomed portion to show details of the 
inner structure. The positions of N and E are indicated, as are the 
anti- solar direction and the heliocentric velocity vector. 
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Fig. 2. continued. P/2010 A2, images from UT 22.3 (f), 23.4 (g), 
25.4 (h) Jan. 2010 using UH 2.2-m, UT 2.3 Feb. 2010 using GN 
(i), UT 19.5 Feb. 2010 using UH 2.2-m G). 



2.1.2. Gemini North Telescope 

The observations were acquired in queue-mode on the 8.1-m 
Gemini North (GN) Telescope on Mauna-Kea, using Gemini 
Multi-Object Spectrograph (GMOS) (iHook et al.ll2004l) . through 



a SDSS / filter (iFukugita et al.lll996h . The detector mosaic is 
composed of three 2kx4k EEV CCDs arranged in a row, result- 
ing in an un- vignetted field of view of 5'. The detectors were 
binned 2x2 resulting in a O'.' 145 projected pixel size. 
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able 1. Elements of P/20 1 A2 and the five known MBCs 
Object a [AU] e i Tj~ 

176P/(1 18401) LINEAR 5:217 OBO 1352 3.172 I inht briclaf= 

133P/(7968) Elst-Pizarro 3.164 0.154 1.370 3.185 " I 

238P/Read (P/2005 Ul) 3.165 0.253 1.266 3.153 .\ 
P/2008 Rl Garradd 2.726 0.342 15.903 3.216 1 

P/2010 R2 La Sagra 3.099 0.154 21.394 3.203 Huc EU5 



P/2010A2 LINEAR 2.291 0.124 5.255 3.588 



2.1 .3. Univertity of Hawai'i 

The observations were performed on the University of Hawai'i 
(UH) 2.2-m te lescope on Maun a-Kea, using the Tektronix 2kx2k 
CCD camera (IWainscoatll20 1 Ol) . The projected pixel size is (7/22, 
and the field of view 7:5. Kron-Cousin B, V, R, I and SDSS / 
filters were used. 

2.2. Data Processing 

The images were corrected for instrumental signature by sub- 
tracting a master bias frame and dividing by a master flat-field 
frame obtained from twilight sky images. Because of the large 
angular size of the object, we have not applied second order, 
"super-flatfield" technique. In the case of the Gemini images, 
the basic processing was performed using the GMOS tasks from 
the Gemini IRAF Package (version 1.9) , which also corrects for 
the geometric distortion and merges the three chips into single 
images. The fr ames were flux calibrated using nightly images of 
lLandoltl([T992h fields obtained over a range of airmasses. 

In order to create deep composite images to assess the ex- 
tent of the dust, off'sets between images were measured and cor- 
rected for using series of field stars. The motion of the comet was 
then compensated for using either the linear ephemeris rates, or a 
full astrometric solution to compute the off'set between expected 
positions of the comet from JPL's Horizon ephemerides. The 
frames were shifted and then combined using either a median 
combination or an average rejecting the pixels deviating the me- 
dian value, in order to remove the stars and background objects 
as well as detector defects and cosmic ray hits. In some cases, 
individual stars were masked and rejected during the combina- 
tion. The resulting 7?-band (or equivalent) images for each run 
are presented in Fig. [21 A zoomed portion of the head of the ob- 
ject is included in each panel in order to show the details of the 
inner structure. 

3. Analysis 

3.1. Description 

Following the nomenclature introduced by iJewitt et aP (1201 Oh 
in their Fig. 1 describing their Hubble observations, the main 
features of the comet are a principal nucleus connected via a very 
faint and narrow light bridge to an arc-shaped dust feature (see 
Fig.O. A second arc-shaped feature crosses the first one almost 
perpendicularly. The main tail extends from these two arcs. Very 
narrow striae along the main tail are emanating from knots in 
the arcs. A low surface brightness fin-shaped structure lies above 
(North) of the main tail. A second low surface brightness more 
diff'use structure (not reported on the Hubble images) extends to 
the SE of the arcs. 

A variety of unsharp masking techniques were used to en- 
hance the tail structures. A filter consisting of a sliding circu- 
lar window of radius n pixels was moved through the 19 Jan. 



NDrLh Uirruiie fe^iurt^ 
\ 

\ 

I 




Fig. 3. Schematic of the main features of P/2010 A2. The back- 
ground image is the 19 Jan. GN image from Fig. Oe. See also 
Fig -S] for enhancement of the features. 



GN composite image and the central pixel in the window was 
replaced by the median of the pixels in the window. The result- 
ing smoothed image was subtracted from the original image to 
remove the background thus enhance the high frequency com- 
ponents. The results for windows of 20, 10 and 5 pixel radii are 
shown in Fig.|4l The 20 pixel radius filtering window brings out 
a broad dust feature at PA ~ 310°, the "North fin", which is 
roughly in the anti- solar direction. The dust modelling described 
below gives some hi nts on the origin of that features, which are 
discussed in Section 13.3.61 Additionally, the shifted individual 
GN images were median combined, and rotated by 8.3° to place 
the tail along a row. 

The rotated image was then shifted by 1 pixel vertically and 
subtracted from the composite to enhance structures along the 
dust tail. This is also shown in Fig. IH panel [d]. Each of the 
condensations of material in the arc are shown to be secondary 
sources of dust. 

The available data-set covers a period of 36 days, during 
which the Earth crossed the object orbital plane (from -2.2° to 
0.6°). The morphology of the comet is remarkably constant over 
that time, implying that the motion of the dust is slow, and that 
the observed dust features have a significant thickness above the 
orbital plane. This, in turn, implies the dust was ejected with a 
non-zero velocity. We searched for evolution in the morphology 
between the various epochs of observation, but given the obser- 
vation interval and resolution of the ground-based images, no 
changes in structure were observed. 

3.2. The nucleus 
3.2.1. Measurements 

The nucleus brightness was measured on the two GN datasets, 
through a series of apertures of increasing radius. The best com- 
promise between enclosed signal and sky background noise is 
found for an 0'.'3 radius aperture. The nucleus magnitude mea- 
sured in that aperture is corrected for the missing flux using the 
star growth profile as a reference. The resulting magnitudes are 
listed in Table [3l 

The radial profile of the nucleus was measured between PA 
-40 and 110° (i.e. the area not aff'ected by the tail). This was 
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Table 2. Observation Log 



Obs. date 






a' 


PsAng"^ 


PsAMV^ 


PlAng^ 


Tel-Ins' 


Who^ 


Obs. 


9 


Seeing 


Comments 


UT 2010 


[AU] 


[AU] 


[deg] 


[deg] 


[deg] 


[deg] 










[arcsec] 




14.2 Jan. 


2.013 


1.040 


5.6 


123.3 


280.1 


-2.2 


NTT 
EF0SC2 


OH 


1200 
600 

2400 
800 


B, 
V, 
R, 


1.6 


Poor seeing 


16.2 Jan. 


2.014 


1.045 


6.7 


117.7 


279.8 


-2.0 


NTT 
EF0SC2 


HH 


900 
600 
1500 
600 


B, 
V 
R, 


1.4 


Poor seeing 


17.2 Jan. 


2.014 


1.048 


7.2 


115.5 


279.7 


-2.0 


NTT 
EF0SC2 


HH 


600 
600 
1200 
600 


B, 
V 
R, 
i 


1.1 




18.2 Jan. 


2.015 


1.051 


7.8 


113.6 


279.6 


-1.9 


NTT 
EF0SC2 


HH 


3600 


R 


0.9-1.6 




19.5 Jan. 


2.015 


1.055 


8.5 


111.5 


279.4 


-1.8 


GN 
GMOS 


OH/GT 


3000 


y' 


0.5-0.7 


Excellent seeing. 


22.3 Jan. 


2.016 


1.066 


10.0 


107.9 


279.1 


-1.5 


UHTelc 
Tek 


BY 


1800 
1750 
7250 
1800 


B 
V 
R, 
I 


0.5-0.8 




23.4 Jan. 


2.017 


1.071 


10.6 


106.8 


279.0 


-1.4 


UHTek 


BY 


3000 


R 


1.0 




25.3 Jan. 


2.018 


1.080 


11.7 


105.2 


278.8 


-1.3 


UHTek 


BY 


4000 


f 


0.7 




02.3 Feb. 


2.022 


1.125 


15.6 


100.4 


278.3 


-0.6 


GN 
GMOS 


OH/GT 


3000 




0.5 




19.5 Feb. 


2.032 


1.262 


22.3 


96.3 


278.3 


0.7 


UHTek 


JP 


1500 


R 


1.2 


Star on object 



Notes: 1, 2: Helio- and geocentric distances; 3: Solar phase angle; 4: Position angle of the extended Sun-Target radius vector; 5: Position angle 
of the negative of the target's heliocentric velocity vector; 6: Angle between observer and target orbital plane, measured from target; 7: Telescope 
and instrument: NTT= New Technology Telescope, GN= Genini North, UH= UH 2.2-m; 8: Observer initials: 0H= O. Hainaut, HH= H. Hsieh, 
GT= G. Trancho, BY= B. Yang, JP= J. Pittichova 9: Total on-target exposure time [sec], and filter. 



Table 3. Nucleus photometry 



UT (exposure start) 


Exp. time 


Airmass 


Filter 


Seeing 


Mag 


M(l,l,Qf) 


M(1,1,0) 






[s] 






n 


r' 


r' 


r' 


2010-01-19 


11:27:03.3 


600. 


1.16 


r 


0.58 


23.93+0.04 






2010-01-19 


11:37:51.3 


600. 


1.19 


r 


0.60 


23.77+0.03 






2010-01-19 


11:48:38.2 


600. 


1.22 


r 


0.63 


23.96+0.04 






2010-01-19 


11:59:25.2 


600. 


1.26 


r 


0.70 


24.01+0.03 






2010-01-19 


12:10:12.1 


600. 


1.30 


r 


0.81 


24.20+0.04 






2010-01-19 


average 






r 






22.34+0.04 


21.74+0.04 


2010-02-02 


06:32:55.1 


600. 


1.09 


r 


0.93 








2010-02-02 


06:43:42.9 


600. 


1.08 


r 


0.93 








2010-02-02 


06:54:30.9 


600. 


1.06 


r 


0.78 


24.03+0.08 






2010-02-02 


07:05:17.8 


600. 


1.05 


r 


0.92 


24.38+0.06 






2010-02-02 


07:16:04.8 


600. 


1.03 
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Note: the nucleus was contaminated by a field star. 



done by converting the (x, y) of each pixel to an angle and a dis- 
tance, then the flux is averaged between the two position angles 
for each distance range. The surface brightness profiles are then 
constructed by iteratively adjusting the asymptotic level of the 
brightness profile so that it reaches zero. The error bars shown in 
Fig. [5] therefore reflect a robust estimate of the sky noise. The av- 
erage profile is compared to that of a field star (measured perpen- 
dicularly to the trailing). The profiles corresponding to the image 
with the best seeing are displayed, normalized to the same peak 
brightness, in Fig. [5j The P/2010 A2 profile shows no signifi- 
cant excess for radii smaller than 3 pixels ((y.'42), then a (noise) 
extended contribution barely above the sky level. 

In order to obtain a conservative estimate of the quantity of 
dust directly surrounding the nucleus, the upper limit on flux ex- 



cess from the nucleus profile with respect to the reference star 
was integrated out to a radius of T\ The upper limit on the flux 
excess in an annulus is defined as the measured surface bright- 
ness of the comet minus that of the comparison star, plus a 1-cr 
error bar, multiplied by the area of the annulus. The limit on the 
flux excess integrated to a radius of corresponds to a mag 
21.1, or < 3% of the flux from the nucleus. The nucleus did 
not present any significant mass-loss at the time of the observa- 
tions. This is further confirmed by the lack of any dust visible 
in the low er-left (SE) quadrant below the nucleus, as discussed 
in Sec tion 13.3.21 This strengthens the conclusion of iJewitt et al.l 
(120101) that there was no cometary activity based on the purely 
geometric variation of nucleus brightness with heliocentric dis- 
tances, withing 0.6 mag, from 25 Jan. until 29 IMay. 2010. 
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Fig. 4. UT 19 Jan. 2010 GN image of P/2010 A2 (Fig. [22), 
enhanced using unsharp masking, with a window of 20 (a), 10 
(b) and 5 pixel radius (c). Panel (d) was produced by shifting the 
image by 1 pixel vertically and subtracting the result from the 
original, after rotating the image to align the tail along the rows. 



The average measured r' magnitudes convert to absolute 
magnitudes /(I, l,a) = 22.34+0.04 and 22.35 ±0.05 for the two 
GN runs. We do not have enough data to estimate the solar phase 
effect. We therefore use the typical S-type asteroid G = 0.15 to 
correct for this effect, assuming the object is either a member of 
the Flora family, or of the most frequent type in that region of 
the solar system. The phase-corrected absolute magnitudes are 
/(1, 1,0) = 21.74 ± 0.04 and 21.55 ± 0.05. The uncertainty is 
however now dominated by the uncertainty of the phase correc- 
tion, which could be off by 0.08 mag if the object is asteroidal but 
not of the S-type (Jedicke et al. 2002), and by 0.10-0.15 if the 
object is a low albedo comet. The magnitude difference between 
the two epoch is likely to dominated by the rotational variability; 
indeed a body this small is likely to have lightcurve amplitude 
of many tenth of magnitude. 

Using the colour equations for the SDSS filters from 
iFukugita et al. l([T99l. the solar magnitude in the Sloan filter is 

= -26.95 (llvezic et al.ll200lh . Assuming a typical S-type as- 




Log Radius [arcsec] 

Fig. 5. Surface brightness profile of the nucleus between PA -40 
and 110°. The individual pixels are represented by small dots. 
The average brightness (between -40 and 110°) in consecutive 
annuli is marked by connected solid disks. The other line is the 
scaled profile of a field star. 



teroid albedo of p = 0. 1 1 , the absolute magnitude corresponds 
then to a radius of 80 and 90 m for the two epochs. The for- 
mal uncertainty on this value is completely dominated by the 
assumption of the albedo and phase correction, and certainly in- 
cludes a variation caused by the rotation of the nuc leus. This is 
in agre ement with the radius of 58-77 m reported by lJewitt et al.l 
(120101) using p = 0.15 (which would convert to 68-90m with 
p = 0. 1 1), and is niarginal ly compatible with the r = 1 10 ± 20 m 
from lMoreno et al.l (l2010h . 

Assuming the body is spherical, with a density of 
3 000 kg m"^ — the a verage of 11 Part henope and 20 Massalia, 
two S-type asteroids (iBritt et al.ll2002h — its escape velocity is 
Ve ~ 0.10-0.12 ms-^ 

If the object were instead a cometary nucleus (with a linear 
phase correction of 0.04 mag/deg and an albedo p = 0.04), then 
the radius would be 120-140 m. In this case, assuming a density 
of 1 000 kg m"^, the escape velocity would be in the range Ve ~ 
0.09-0.11 m s"^. As discussed later, we favour the asteroidal 
nature of the object. 

Photometric measurements in the other filters have insuffi- 
cient signal-to-noise ratio to produce meaningful colours. 

3.2.2. Thermal model 

In order to examine the possibility of ice sublimation being the 
driver of the observed dust ejection features, we treat the ob- 
ject as a comet nucleus. As such, the question is how long will 
ice survive in the interior of such a small object at such a small 
heliocentric distance. For the purpose of modelling the thermal 
evolution of the nucleus, we assume a similar 1-D m odel and 
numerical i mplementat i on to t hat described in detail by lPrialni^ 
(199^ and lSarid etaP (l2005h. Further disc ussion of the input 
physics can be found in lPrialnik et al.l (|2004|) . 

A comet nucleus is generally portrayed as a porous aggregate 
of ices and solids (IWeidenschillingI 120041) . The internal struc- 
ture can be modeled as an agglomeration of grains made of 
water ice (and perhaps other minor volatile compounds) and 
solids (silicates, minerals) at some mixing ratio, with a wide 
spread in the size distribution of the components. The ice more 
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volatile than water are neglected in our treatment as any such 
ice becomes depleted during a long term evolution, even down 
to the core of much larger objects (see Prialnik & Rosenberg 
I2OO9I) . Furthermore, if the object originated in the asteroid belt, 
it is more likely that H2O ice would be crystalline and that 
the initial abundance of ni ore volatile ices would be minute 
(e.g. lMeech & Svorenll2004l) . We solve the set of coupled, time- 
dependent equations of mass cons ervation (dust, ice a nd vapor 
components) and heat transfer (see lPrialnik et al.|[2004 for a de- 
tailed description). The boundary conditions for these evolution 
equations are vanishing heat and mass fluxes at the center, van- 
ishing gas pressures at the surface and a requirement for energy 
balance at the surface: 

F(R) = e(TT(Rjf + 

P,,,(T) ^^nRJH - (1 - A)^^ cos^, (1) 

where e and A are the emissivity and surface albedo of the nu- 
cleus, Lq is the solar luminosity, cr is the Stefan-Boltzmann con- 
stant, Rg the molar gas constants and yu, Pyap and H are the molar 
mass, vapor pressure and sublimation latent heat for water. This 
condition on the surface flux, F(R), depends on the Temperature 
T , heliocentric distance r and local solar zenith angle ^ and is 
calculated for each time step. The first term corresponds to the 
re-radiated heat, the second to the sublimation, and the third the 
absorbed solar radiation. Without any additional information, we 
assume a fast rotator with the sub-solar point illuminated at per- 
ihelion. 

We consider two cases for this object: "Case I", which is 
similar to an S-type asteroid, and "Case 11", which is taken as 
a comet nuclei. Using the assumptions and measurements de- 
scribed in the previous section. Case I has a radius, bulk density 
and albedo of 85 m, 3000 kg m"^ and 0.1 1, respectively. Case II 
has a radius, bulk density and albedo of 130 m, 1000 kg m"^ 
and 0.04, respectively. If we assume the solid component (both 
the nucleus and the dust) to have specific density similar to that 
found for the d ust particles of comet Wild 2, ~ 3400 kg m"^ 
(iKearslev et al.1 12009), we can constrain the initial ice abun- 
dances for the thermal evolution calculations (see lPrialnik et al.l 
[2QQ8h . These values are taken as 0.05 and 0.5 (with correspond- 
ing initial porosities of 0.01 and 0.3) for Case I and Case II, re- 
specti yely. All other physical param eters are taken to be similar 
to the lPrialnik & Rosenberg (l2009l) . 

Our simulations of the nucleus evolution are run until the wa- 
ter ice component is depleted from the interior. For Case I ("as- 
teroid"), there is no more water ice present after ~ 3 x 10^ yr, 
while for Case II ("comet") water ice becomes negligible and is 
buried deep under the surface after several 10^ yr. This is shown 
in Fig. [6^, where the sub- surface water ice depth is plotted as a 
function of time. Beyond this point in time, the water ice com- 
ponent is rapidly depleted, as the innermost layers reach a tem- 
perature of 200 K and 1 80 K, for Case I and Case II respectively. 
This is shown in Fig.[6j3, where the evolution of the central tem- 
perature is plotted, for both cases, up to 10^ yr. The penetration 
of heat into the interior and the rise of the innermost temperature 
is slower for Case II and the final temperature reached is lower 
throughout most of the evolution time. Case I heats up faster be- 
cause it is less porous (bulk density is higher) and has a much 
larger mass fraction of dust (silicates and minerals), so heat dif- 
fusion is less attenuated by sublimation. 

In summary, considering a comet-like and an asteroid-like 
cases for the composition and structure of the nucleus, our ther- 




10^ 10"" 10^ 10^ 10^ 

Time [yr] 




10^ 10' 10' 10^ 10^ 
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Fig. 6. Evolution of the water ice sublimation front and the cen- 
tral temperature for the nucleus of P/2010 A2. The x-axis rep- 
resents the time [yr] in log scale since the placement of the ob- 
ject in its current orbit, the y-axes represent: a. The sub-surface 
depth where crystalline ice survives; b. The central temperature. 
Case I ("asteroid") is plotted in black, while Case II ("comet") 
is plotted in gray. 

mal modelling eff'orts indicate that the interior reaches a temper- 
ature over 180K in a few 10^ to a few 10^ years. The nucleus 
would therefore be fully depleted of water ice — and any other 
more volatile ice — in that time scale. The orbit of P/2010 A2 
is typical of an inner Main Belt asteroid, possibly member of 
the Flora family, with no indication that it would have been very 
recently injected or captured from a location more distant from 
the Sun. The survival depth of water ice becomes greater than 
~ 10% of the object's radius in less than a million years, al- 
most independent of the initial choice of model parameters. This 
model therefore indicates that ice sublimation could not play any 
role in lifting dust from the surface of the object or dragging any 
substantial amount of grains from the interior. In what follows, 
we will therefore assume that the nucleus of P/2010 A2 has the 
characteristics of a Flora family, S-type asteroid, as detailed in 
this section and the previous one. 

3.3. Dust and the tail 

3.3.1. Dust dynamical modelling 

The dust tail was mo deled using the Finson-P robstein (FP) 
dust dynamical method (iFinson & ProbsteinllT968b . modified by 
iFarnhamI (11996b . 

Dust is extracted from the nucleus, by gas drag in traditional 
sublimation-driven cometary activity, or via direct ejection in the 
case of an impact. The dust is decoupled from the gas flow (if 
any) within a few nuclear radii, and is no longer influenced by 
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the nucleus gravity, which becomes neghgible. The motion of 
a refractory dust grain will then be defined by its initial veloc- 
ity and aff'ected by the solar gravity and radiation pressure, both 
acting in opposite radial directions. The net result of the two in- 
fluences may be thought of as a reduced gravitational force act- 
ing upon the dust. The parameter p is defined as the ratio of the 
radiation pressure force to the gravitational force, and is given 
by 

yS = 5.740 X 10-4^ (2) 
pa 

for grain of radius a [m] and density, p=3 000 kg m"^. Q^^ is 
the radiation pressure efficiency (typically 1-2, depending on the 
material scattering properties). 

The FP method involves computing the trajectories of dust 
grains ejected from the nucleus at velocity v and under the in- 
fluence of solar gravity and radiation pressure. A collection of 
particles of different p ejected at a particular time follow trajec- 
tories called synchrones, while a single particle size {J3) ejected 
over a range of times follow trajectories called syndynes. The 
combined set of curves maps out the expected ejected dust dis- 
tribution. The extent to which the actual dust surface brightness 
correlates with the syn-curves may be used to infer which grain 
sizes are present, the onset an duration of activity, and any initial 
ejection velocity. 

In order to investigate the origin of the dust tail, rather than 
starting with an entirely generic model which would have a large 
number of loosely constrained free parameters, we proceeded 
by steps, focusing on some aspects of the tail and making some 
simplifying hypotheses. 

3.3.2. Tail morphology and duration of the activity 

Overlaying the syndyne and synchrone curves onto the image of 
the comet (see Fig. (T]), it is clear that the region corresponding 
to dust emission during the 30 days before the observations is 
completely empty of dust. Together with the negligible contri- 
bution from the dust to the near-nucleus area, this confirms that 
any cometary activity must have stopped at least several weeks 
before the observations. 

The light bridge connecting the nucleus to the tail matches 
the po sition expected for very large dust grains. iJewitt et aP 
(l201Qh measured the orientation of the tail on the high-resolution 
HST images, concluding that there was a short period of dust 
emission between Feb. and Mar. 2009. 

However, while the orientation of the tail corresponds well to 
the direction of that synchrone, its position is offset from what 
is expected from a dust emission a zero speed from the nucleus. 
Also, the position of the northern fin is not matched by any syn- 
dyne, and the constant appearance of the tail while the Earth 
passed through the orbital plane indicates a more complex, 3D 
geometry and dust emission with non-zero velocities. 

We will now further investigate the dust release character- 
istics. Adding up the scattered light from a distribution of dust 
particles emitted from the nucleus over a range of times as it 
moves away from the object, one can build synthetic images of 
the corresponding tail. Comparing the surface brightness of the 
synthetic tail to the real data, observational data can be inverted 
and information about the grain size distribution, ejection veloc- 
ity and onset and cessation of activity can be obtained. 

Various general dust emission models were explored to in- 
vestigate the overall morphology of the tail. The first model in- 
cluded emission for a long duration of time, about a year from 



the time of observation to three months before the observation. 
The second model included only one month of emission begin- 
ning a year prior to the observation. Finally, a third model only 
included a single burst of emission one year before the obser- 
vation. The emission takes place at the sub-solar point, with a 
constant dust production during the considered period. The dust 
size range and distribution and emission velocities were adjusted 
so that the reconstructed tail reproduces as well as possible the 
observations. The parameters are adjusted to reproduce the ori- 
entation (measured by its position angle) and width (measured 
perpendicular to the tail) of the tail, and the surface-brightness 
distribution along the tail. The actual position offset between the 
tail and the nucleus is not measured at this stage; it is considered 
as an effect of the global/average emission velocity, and will be 
dealt with later. 

The best parameters for each model are listed in Table IH the 
corresponding reconstructed images are shown in Fig. [8] and the 
surface brightness profiles are displayed in Fig.|9l These profiles 
were obtained by averaging 52 pixels perpendicularly to the tail 
on an image rotated by 8° so to have the tail along the X direc- 
tion. The first two models failed to reproduce the overall width, 
shape and/or brightness profile of the tail, even for the best ad- 
justment of the parameters. While it may not be obvious from 
Fig. [HI the light distribution along the length of the tail is very 
different for models 1 and 2 compared to the reference image, as 
seen in Fig.|9l 

The third model, with a single burst of emission, reproduces 
well these shape, width and profile of the tail with the particle 
sizes in the a = 1.5 to 20 mm range, with a size distribution 
proportional to a~^-^. The low end of the range is set by the edge 
of the image, and the upper end corresponds to particle with a 
minimal sensitivity to the radiation pressure. 

The parameter values for Model 3 are tightly constrained in 
the adjustment process: power laws shallower (-3.75) or steeper 
(-3.25) fail to reproduce the profile of the tail. The best range 
of speeds of the particles was 0.20 to 0.30 m s"^ with smaller 
particles having higher velocities. The range of particle speeds 
that still develop an acceptable model is only 0.10-0.40 m s"^ 
The total emission can vary by less than a factor of 2 and still 
acceptably reproduce the observed surface brightness. The total 
amount of mass ejected (assuming a density p = 3 000 kg m"^ 
and albedo p = OAl) was about 5x10^ kg, again within a factor 
of 2. 

In summary, reproducing the overall morphology of the tail 
suggests that the dust was emitted during a single burst about one 
year before the observations. The dust grain distribution follows 
a power law in -3.5 in the 1.5-20mm range, emitted with slow 
velocities in the 0.20-0.30 m s"^ range. 

We al so note tha t , for a water-ice sublimation-driven dust 
emission, iDelsemm^ (Il982h obtains for an object at r ~ 2 AU 
ejection velocities in the 0.3-0.4 km s"^ range, i.e. over 3 or- 
ders of magnitude larger than the velocities obtained in the FP 
modelling. However, the assumptions made by Delsemme are 
not necessarily valid her e, so we will not use this argument. 

iMoreno et al.l (l2010l) analyzed a set of images obtained be- 
tween 14 and 23 Jan. 2010 using a different modelling technique, 
assuming particles ranging between 0.01 and 10 mm with a den- 
sity of 1000 kg m"^ and an albedo of 0.04. Their best fit suggests 
dust emission lasted over several months, ending just after peri- 
helion (Dec. 2009), with a complex emission function over that 
period. While their model was richer than ours in terms of the 
complexity of the emission pattern, our first FP analysis cannot 
support that there was dust emission that late (see Fig. (T]). Our 
model that is closest to their best fit is Model 1 (although ours 
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Fig. 7. Syndynes (labeled with (3 = 0.01-0.05) and synchrones (labeled in days since emission, 18-345) for dust grains emitted from 
the nucleus with zero velocity. The background image if from Figl2le. 



Table 4. Three FP models 



Model 1 2 3 

Emission period, 

starting on 10 Feb. 2009 250 days 30 days 1 day 

Particle size range (mm), a 0.3-3 1.5-100 1.5-20 

Particle size distribution power -1 -3.9 -3.5 

Emission velocities (m/s) 0.20-0.30 0.10-0.20 0.20-0.30 

Fidelity of the model Poor Very poor OK 



mating grains in traditional comets. It is also a warning that ba- 
sic FP modelling should be applied with extreme care. In the 
case of P/2010 A2, however, the situation is actually simpler. 
Considering the bright knots in the arcs as the secondary dust 
grain source giving birth to the corresponding striae, the straight 
striae emanating from the knots follow synchrone curves. The 




100 200 

Position [pixels] 



Fig. 9. Photometric profile of the reference image (bold line) 
and the three models. Model 1 was shifted vertically for clarity. 



range of emission epoch compatible with the observations is 
measured esitimating (by eye) the range of position angle that 
fit that of the striae. This range corresponds to emission taking 
place 340^20 ^^y^ before the observations, i.e. in agreement with 
the epoch of the main event. We conclude that each stria corre- 
sponds to dust grains with a distribution of JS, ejected from the 
nucleus at the same time of the main dust release, with the same 
velocity as the material in the knot leading the stria. The sur- 
face brightness distribution as a function of the distance to the 
arc constitutes therefore a "dust grain radius spectrum" of the 
material. 



Note: the emission duration is set, the other parameters are adjusted to 
reproduced the observed tail. 

uses a constant emission). Our Mode l 3, with an burst emis- 
sion, reproduces the tail much better. iJewitt et aP (1201 Ol) con- 
cluded from the evolution of the tail position angle over 25 Jan.- 
29 May. 2010 that the dust had been emitted in a brief burst 
around 2 Mar.2009, plus or minus a few weeks. [Snod grass et al.l 
(I2OIOI) , taking advantage of the diff'erent viewing geometry from 
the Rosetta spacecraft, concur that the emission was very short, 
further constraining the emission date to 10 Feb. 2009, plus or 
minus 5 days. 

Model 3 nicely reproduces the overall morphology of the 
tail, in particular its orientation and width, and the detach- 
ment from the nucleus and the overall surface brightness pro- 
file. However, it fails to reproduce the sharp edges and features 
visible in the head. Using the results from this model, we will 
further investigate the hypothesis of a single burst of emission. 

3.3.3. Quantity of dust in the tail 

The presence of striae in the tail is usually a sign that sec- 
ondary dust emission is taking place, for instance from subli- 
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Fig. 8. Synthetic image of P/2010 A2, a. for 
the best fit with a long period of dust emis- 
sion (Model 1, in Table |4]). b. for the best fit 
with a month-long period of dust emission 
(Model 2). c. for the best fit with single burst 
of dust emission (Model 3). d. Reference im- 
age, same scales. This is the same image as 
Fig. 2.e. 



b. 




d. 



We wish to estimate that dust size distribution, as well as 
the total quantity of dust in the tail, using the sky- subtracted, 
flux calibrated composite image from 19 Jan. 2010 (the deepest 
image). Eor each pixel we want to know which fraction is cov- 
ered by dust (filling factor /), multiplied by the albedo (A). Af 
is a convenient dimension-les s equivalent to the su rface bright- 
ness. Using the formalism of lA'HearnetaD (11984 we convert 
the measured flux in A/, ie. 



Af- 



1 2AkmrAU 
\ Pkm / 



(3) 



Akm is the geocentric distance in km, tau the heliocentric dis- 
tance in AU, Pkm the linear size of the aperture at the distance of 
the comet (in km), Fo^ the measured flux of the comet, and Fq 
the flux of the Sun at lAU. Converting Eq.[3]for square pixels {p 
arcsec), and expressing Fo</Fq in CCD ADUs, it becomes 



Anr- 



(4.848 X 10-6;?)^ 



^Q-0.4(ZP-/:(z-l)-Mo)^j3U^^ 



(4) 



where ZP is the zero point, z the airmass and k the extinction 
coeflftcient (in mag/airmass; Gemini's ZP is expressed for 1 air- 
mass, ZP = 28.20 from Section [T2b . 

In order to count the number of dust particles in each pixel, 
we estimated the size of the particles, assuming that all the dust 



in the tail was emitted during a single burst that took p lace 
on 10 Eeb. 2009 (using the date of ISnodgrass et al.l (l2010li . see 
13.3.21) , with a zero velocity from the large chunks of material 
that ended up located in the arcs. The apparent distance from 
the arcs to the particle is then a measurement of J3, the conver- 
sion factor is estimated running our EP code for particles over 
a range of values of jS for an emission date in Eeb. 2009. Using 
the relation between B and the characteristics of the particle from 
lEinson & Probsteiiil (Il968l) . we convert in the radius of t he par - 
ticles, a, using Equation [2l We use Qpr = 1, as iMorenol (l2009l) . 
for large absorbing grains (iBurns et al J[T97l . We assume a den- 
sity of 3 000 kg m"^ (see Section [3". 2. lb . Eigure [TOl displavs the 
corresponding sizes, ranging from few tenths of millimetres at 
the edge of the fr ame to few mill i metres close to the arcs. Using 
the formalism of I Agarwal et al ] (I2007h . the dust size distribu- 
tion inferred from the surface brightness profile along the tail 
is a power law in the a = 1-20 mm range, with a number of 
particles proportional to a" (dififerential size distribution), with 
n = -3.44 ± 0.08. This value is in agreement wit h that modeled 
in the previous section, and with tho s e repo rted by lMoreno et al.l 
(201Qb, n = -3. 4 ± 0.3. lJewitt et aP (l2010h . n = -3.3 ± 0.2 and 
ISnod grass et al. (l20Tol) . n = -3.5. A value close to -3.5 is sig- 
nificant because this is in agreement with the empirical value 
for the distribution of particle sizes resulting from a collision. 
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and for a collisionally relaxed population (iDohnanyill 19691) . The 
value for the dust in the tail of P/2010 A2 is slightly steeper than 
the size distribution measured with Hayabusa for the boulders 
larger than 5 m on the surface of (25143) Itokawa, n = -3.1 
(iMichikami et al.ll2008h . 

The total Af is measured in a series of thin rectangular aper- 
tures as long as the width of the tail and covering it from the arcs 
until the edge of the frame. The total Af is converted into a num- 
ber of particles using the radius a of the particle at that distance 
from the arcs as computed above. 

The total mass and volume of dust is integrated by adding 
the contribution of each rectangular aperture starting from the 
arcs. The mass integrated along the tail, as well as the volume of 
the corresponding sphere are displayed in Fig. [TOl While there is 
still a significant amount of dust in the tail extending beyond the 
edge of the frame, the asymptotic behaviour of the growth curve 
indicates that a very large fraction of the mass is present in the 
visible part of the tail. With these assumptions, we estimate that 
the tail contains at least 8 x lO^kg of ejected dust, which could 
be re-assembled into a sphere of radius ~ 40 m. 

Jewitt e t'aP (120101) reached similar estimates based on the 
analysis of their Hubble images: 6-60x10^ kg, corresponding 
to a sphere with a radius of at least 17-36 m. The latter corre- 
sponds to a distance of 20'' from the arcs on Fig. \W[ which is 
about where the S/N starts to degrade on their Hubble images: 
we consider their larger estimate in agreement with ours. The 
total ejected mass is likely to be higher, as these values poorly 
estimate the size and number of large particles in X-shaped arcs. 

3.3.4. Envelope of the tail 

A dust emission model was developed to investigate the possible 
cause of the position of the tail behind the X-shaped arcs, and 
its distinct detached geometry. For this investigation, rather than 
trying to reproduce the light distribution in detail using a general 
fit of the dust velocity vector Vd and JS distributions, a greatly 
simplified model was implemented as follows. The details of the 
model, of its implementation and of its results are beyond the 
scope of this paper, and constitute a separate, companion paper 
(iKlevna et aLlboilh . 

The dust emission is r epresented by an insta ntaneous burst 
on 10 Feb. 2009 (date from lSnodgrass et al.l201Ql see lTT^ . The 
first simplification represents the emission geometry with a cone 
(defined by its half-opening angle ac and axis orientation), on 
which dust grains characterized by a distribution of j3 are emitted 
with a distribution of velocity vj. The motion of the dust grains is 
computed taking into account the Sun's gravity and the radiation 
pressure to determine their final positions in the image plane at 
the time of the observations. 

The second simplification is to only consider whether the 
dust grains fill the envelope of the observed tail, completely ne- 
glecting the light distribution within that envelope. The envelope 
is defined visually on the images of 19 Jan. (Figl^e), so as to 
enclose the main tail. Several diff'erent plausible envelopes, esti- 
mated by eye, were defined to test the dependency of the result 
on the choice of the envelope. This allows us to neglect the ac- 
tual distributions of J3 and va: only the minimum and maximum 
values are relevant to the model. The natural minimum veloc- 
ity and minimum dust size are zero; the dust size is assumed to 
be independent of the velocity (this is supported by experimen- 
tal studies of impact onto porous, regolith-like material, which 
are appropriate for gravity- or weak cohesive strength-controlled 
craters which, as we discuss later, seems to be the case here). 
Finally,the only diff'erence between a filled emission cone and a 



hollow one is the presence of a cap at the open end of the cone, 
we can therefore further simplify the model considering a hollow 
cone. 

This leaves only four free parameters: the two coordinates 
defining the orientation of the cone axis, the cone half-opening 
a, and the maximum dust velocity Vd,max. The merit function of a 
set of parameters is computed by assigning penalty to dust grains 
ending up outside of the envelope of the tail, and to pix els within 
the tail that contain no dust grain. iKlevna et al.l (l201ll) describes 
the exploration of the four parameter space, the method used and 
the sensitivity to the definition of the envelope. Only a fairly nar- 
row subset of the parameter space can reproduce the observed 
tail envelope. The best model is a cone (either full or hollow) 
with a half-opening of 40°, pointing below the orbital plane and 
on the forward side of the comet, with a maximum dust veloc- 
ity of 0.2 m s"^. The direction of the emission is strongly con- 
strained by the position of the tail with respect the nucleus and 
by the sharp edge of the head. 

3.3.5. Arc-shaped features 

In order to investigate the origin of the arc-shaped features (see 
Fig. [3]), ano ther model was imp lemented; it is also described 
in detail in iKleyna et'aP (1201 lb . Each arc appears as a one- 
dimensional feature; they are not parallel to the synchrones 
for the date of the dust emission (see Fig. [71 position angle 
near 278°). They can therefore not be caused by the radiation 
pressure-driven spread of particles with a range of jS; they must 
therefore correspond to very large particles very weakly sensi- 
tive to radiation pressure (i.e. with J3 = 0). 

A large number of particles wer e modeled, with an em is- 
sion date on 10 Feb. 2009 (date from lSnodgrass et aDl2010l see 
13.3.2b , sampling the emission direction and velocity space. The 
motion of these particles was then computed until the date of the 
observations, and their position in the image plane computed. 
Those landing on position of the arcs are marked. Only particles 
with a ejection velocity Vd > 0.20 m s"^ can contribute to the 
X-shaped arcs, but this model cannot set an upper limit. Indeed, 
dust grains with higher emission velocity can be emitted and 
end up further away from the nucleus at the time of the observa- 
tion, but still appearing in the same position in the image due to 
foreshortening. Test particles from two regions in the (velocity, 
direction) space end up in the position of the arcs. The arcs could 
therefore have been produced by dust grains from one of these 
region, the other, or a combination of both. 

However, one of these regions is found to match a section 
of the surface of the cone modeled to describe the overall tail 
in the previous section: the direction of emission and velocity 
are identical. This suggests that the same hollow cone (with its 
orientation, opening and emission velocity) describes both the 
overall tail geometry as well as the arcs in the head of the tail. 
The former corresponding to small dust grains subjected to radi- 
ation pressure, the latter corresponding to a fraction of the large 
dust grains populating the propagated cone. Only a fraction of 
the large dust grains are visible as arcs, as the geometry at the 
time of the observations gives a foreshortened side view of a 
segment of the cone. The other large dust grains appear diluted 
over the tail. 

Observations from another point of view should show other 
features fortuitously aligned alo ng that other line of sight. 
Unfortunately, the Rosetta images (ISnodgrass et al do not 

have a sufficiently high resolution to reveal them: all the region 
presenting structure was contained in one single pixel of the 
Osiris Narrow Angle Camera on Rosetta. In this context, the 
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Fig. 10. The labeled curves show, as a function of the distance 
to the arc-shaped features in the head of the tail, the radius of 
the dust particles [mm], the integrated mass in the tail out to that 
distance [kg] and the radius of the sphere of equivalent volume 
[m] . Refer to the text for the assumptions that were used to derive 
these values. 



striae visible in the tail (see Fig. [3] and |4]) simply correspond to 
irregularities of the emission along the cone. 



3.3.6. Fin-shaped features 

The position of the North fin, above and beyond the asymptotic 
direction of the synchrone in Fig. [71 suggest an out-of-plane fea- 
ture, in turn indicating an out-of-plane ejection velocity vector. 
It can be explained as a consequence of the ejection cone de- 
scribed in the previous section. The cone itself defines the left 
edge of the fin, the smaller dust being driven back to create the 
body of the feature. The SE fin could not be reproduced. 



3.4. Alternative emission processes 

When evaluating the characteristics of the nucleus, in Section 
13.2.21 we concluded that no water ice or more volatile species 
would have survived for more than a few million years in the 
body on its current orbit. This makes cometary activity a very 
unlikely cause for the observed dust emission, as it would require 
a mechanism to inject the object on its current Flora-like orbit in 
the recent past. 

The FP m odelling of the tail (Section [3. 3. lb and other stud- 
ies (iJewitt et al. 2010; Snodgrass et al. 2010) firmly indicate the 
emission duration was very short, in the form of a burst that took 
place around Feb. 2009. The more detailed dynamical dust mod- 
els presented in the previous section suggest this emission took 
place on a cone with an half-opening a = 40°. As we will dis- 
cuss in the following section, this strongly suggests that the dust 
emission corresponds to an ejecta cone resulting from an impact 
with a smaller body. In this section, we shall briefly consider 
other alternative dust emission processes: rotational spin-up and 
electrostatic dust levitation. 



Rotational spin-up : SmaU objects in the main belt can 
be driven towar d rotational instability by radiation torque 
iRubincamI (l2000l) . Radiation torque can spin-up or spin-down 



the object depending on its shape and the nature of the ob- 
ject, and the numerous small asteroids that are either fast 
and slow rotator s sugg est this efl'ect is indeed acting widely 
(iPravec & Harrisll2000l) . However, simulations indicate that the 
ejection of material from a spu n-up object occurs along the equa- 
torial disk (IWalsh et al.ll2008l) . i.e. with a geometry that is very 
difl'erent from what is observed here. These numerical simula- 
tions were done for particles larger than dust. However, the ge- 
ometry of the large grains observed in this object seems com- 
pletely incompatible with the rotational spin-up lifting. 



Electrostatic lift: Spacecraft have observed fine particles levi- 
tated at velocities near 1 m s"^ at the Moon's terminator caused 
by the buildup of a potential difi'erence between illuminated and 
shadowed regions (Ide Bibhas & Criswellll 19771) . On small aster- 
oids, this eff'ect can launch dust above the local escape velocity 
(Lee 1996), creating visible mass loss. However, there are thou- 
sands of asteroids similar in size to the MBCs that do not eject 
observable dust. Furthermore, the FP analysis of the dust around 
P/2010 A2 indicates it is constituted of large dust grains, which 
could not be lifted up by electrostatic force. 

3.5. Physical implications of ejecta geometry 

One remaining source of the observed release of dust is the exca- 
vation of material from an impact on the nucleus. In this section, 
we demonstrate to first order that an impact event is consistent 
with, and could explain, the three main requirements from obser- 
vations: the duration of release, velocity distribution and prove- 
nance, and amount of mass. We restrict our analysis to a canon- 
ical vertical impact here, although impacts at normal incidence 
rarely occur on planetary objects; instead, the vast maj ority oc- 
cur obliquely for both sphe rical ( Gilbert 1893; Shoema keill963h 
and asteroidal-type bodies (ICheng & Barnouin-Jhalll999l) . 

In a hypervelocity impact event, strong Shockwaves and sub- 
sequent rarefaction waves set material in motion, excavating 
and ejecting a portion of the displaced target up and out of the 
crater; the remainder of the displaced mass remains compressed 
in the target. In general, the velocities of main-stage ejecta can 
be described by a power-law relationship with either time or 
launc h position, as dimensionally predicted (e.g. jHousen et al.l 
'198 3|), and demonstrated by a number o f expe ri mental (e.g., 
Mc Getchin et al.l 119731: ISchultz & GauE Il979l: iPiekutowskil 



198 01: lOberbe ck & Morris on! 119761: iPiekutowski et al.l Il977[ 



ICintala et al B999; An derson et al. and numerical studies 



(e.g. lWada et a l. 2006). The ejection angles of ejecta are a func- 
tion of the material properties of the target. In granular target ma- 
terials (like regolith, as expected on an asteroidal surface), ejec- 
tion angles (referenced from local horizon tal) (A = 90° -de range 
fro m - 30° to - 5 0° (Cintala et al."" 1999'; "Anderson et alj|200|; 
An derson & Schult z 2005; Hermalvn & Schultz 2010) with an 
average ac ~ 45° (IHousen et al.1 119831) and evolve through- 
out crater growth. In cohesive targets (such as solid basalt or a 
welded materia l), ejec tion angles tend to be considerably higher 
(i// ~ 70° lGaultiri973h . Additionally, highly porous material ex- 
hibits a high-ejection angle component ( Schultz et al. 2007). Our 
finding of an ac = 40° hollow cone is in agreement with an 
impact into an unconsolidated, granular target. The progression 
from lower ejection angles (larger ac) for faster velocity material 
to higher ejection angles (smaller ac) at later times and slower 
velocities is consistent with the ejection angle evolution in ex- 
peri ments and computations for unconsolidated granular m ate- 
rial (iHermalvn & Schultzll2010l:lHermalvn & Schultzll2011l) . 
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Crater growth is arrested either when the ejecta velocities 
are insufficient to launch material above the crater rim (gravity- 
controlled) or when the shock strength drops below the yield 
strength of the material (strength-controlled). The cutoff in ve- 
locity solution space of ~ 0.10 m s"Ms coincident with the es- 
cape velocity of the body. Therefore, any material ejected under 
this velocity component could have returned to the surface of 
P/2010 A2 and no record of this material would be apparent in 
the images. This behavior would be consistent with a gravity- 
scaled crater. Alternatively, if the surface exhibits a degree of 
material strength, crater growth could be arrested with virtually 
no ejecta launched below the 0.10 m s"^ cutoff. We note that 
although it is impossible to separate the effects of gravity and 
strength controlled crater arrestment without watching the evo- 
lution of the ejecta curtain. Since both gravity and strength are 
inferred to be weak, both probably play a role in the termination 
of growth. It is possible to define a theoretical gravity-controlled 
crater radius as a limiting size since ejection velocities during 
main- stage growth will be unaffected by the difference in termi- 
nation processes. 

In Section [3". 3. 31 we concluded that the visible ejecta corre- 
sponds to a sphere with a radius ~ 40 m. iJewitt et aP (1201 Ol) 
estimated in the 17-36 m range. In both studies, the dust par- 
ticles were found in the 2-20 mm range, which further supports 
the excavation of fine regolith-like material in a gravity- scaling 
regime. To approach the size of the crater, we look to the ejected 
mass as a constraint. The total amount of mass ejected from 
a crater tends to only be approximately 50% of the total dis- 
placed mass fe.g. lStoeffler et al.lll975l:ICroftlll980l) . The cumula- 
tive (and incremental) ejected mass above a certain velocity sim- 
ilarly follows a power-law relationship with velocity such that 
the vast majority of material is ejected at low velocities near the 
crater rim. The form of this equation fo r purely gravity- scaled 
craters is (as given in lHousen et al ] fT983h : 



' R= 25 m 

- R= 50 m 

- R= 100 m 
R= 150 m 

- Measured Volume 



\^^8R 



-1.22±0.02 



(5) 



where Vg is the cumulative volume of ejecta with higher ve- 
locity, R is the final gravity- scaled crater radius, v is the veloc- 
ity of ejecta, and g is the gravity on the surface; the exponent 
and coefficient are fit to data from small-scale explosion and 
impact experiments into granular targets. Using the higher den- 
sity values defined above, a crater approximately 100m in radius 
is required to eject sufficient material above escape velocity to 
match the observed minimum ejected mass (see Fig. [TT1) . The 
total time to eject the observed material would be on the order of 
seconds. A strength-co ntrolled crater follows a similar relation- 
ship (see iHousen et al ] ([T983i) for details), but would require a 
much larger crater to eject the measured volume of material for 
anything other than extremely small values of material strength 
Y. Considerably larger craters are probably not physically possi- 
ble on such a small body since even the gravity- scaled limiting 
radius is on the order of the radius of the body itself. However, 
there are other bodies that exhibit craters or large-scale depres- 
sion features of radii approaching that of the body itself; e.g., 
Mathilde, Phobos and Deimos, Gaspra, and even the Moon. 

One possible mechanism that could help explain the develop- 
ment of such large craters without catastrophically disrupting the 
body is an oblique impact, which decreases the peak sho ck pres- 
sures while still e xcavating a large amount of material (Schultz' 
1997b ICheng & B arnouin-Jha 1999; Schultz, P H. & Crawford, 
20111) . This process could also explain the enhancement of 
higher velocity material (> 0.25 m s"^) in the images, since 




velocity (m / s) 

Fig. 11. Cumulative volume of ejecta with higher velocity for a 
set of crater sizes in the gravity- scaled regime. Measured vol- 
ume of material is indicated by solid black line. A crater approx- 
imately 100m in radius is required to excavate sufficient material 
above the escape velocity of the body to match observations. 

ejecta downrange (i.e., in the direction of impact) is en- 
hanced in mass and ve locity throughout much of crater growth 
(I Anderson et al. I [2003b and may appear strongly foreshortened 
in the images due to perspective. The body may also be fairly 
porous, which will serve to further attenuate the peak shock pres- 
sures. 

4. Summary and discussion 

The dust tail of P/2010 A2 was observed and analyzed so to cast 
some light on the possible cometary nature of that object, as its 
appearance and orbit could make it a member of the small family 
of Main Belt Comets. The main results and conclusions of this 
study are listed hereafter. 

- Observations were acquired with GN, NTT, UH 2.2-m from 
14 Jan. until 19 Feb. 2010. The tail did not show any measur- 
able evolution over that time span, which included the cross- 
ing of the orbital plane, which indicates the tail is an out-of- 
plane 3D structure, suggesting the dust was released with a 
significant orthogonal velocity. 

- The nucleus, detached from the tail, appears not to be sur- 
rounded by dust (at most 3% of the light corresponds to 
near-nucleus dust); dynamical dust models also indicate that 
no recent dust is present, suggesting no cometary activity at 
the time of the observations. The absolute magnitude of the 
comet was estimated to SDSS /(1, 1, 0) = 21.74 ± 0.04 on 
19 Jan. 2010 and 21.55 ± 0.05 on 2 Feb. Assuming an albedo 
p = 0.11, a value typical for S-type asteroids, this converts 
into a nucleus of radius r = 80-90 m. Using a density of 
p = 3 000 kg m"^ (also typical of a S-type asteroid), the es- 
cape velocity of the body is = 0.10-0.12 m s"^. 

- A thermal model of the nucleus indicates that water ice (as 
well as any more volatile ice) would not survive more than 
a few million years in the object on its current orbit. As we 
have no reason to suspect a recent insertion on the current 
orbit, this rules out ice sublimation as the source of the ob- 
served dust. Rotational spin up and electrostatic lifting were 
also rejected as possible source for the dust tail. 
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- A Finston-Probstein-type dynamical dust modelling of the 
morphology, light distribution and detachment of the tail in- 
dicates that a short burst of dust emission (a day long or 
less) represents best the observed tail, with an ejection ve- 
locity in the range = 0.20-0.30 m s"^. The ejection date 
is in a greement with ot her studies of this object (I Jewitt et aP 
I20T Ol: ISnodgrass et al.ll2010l) : we continue this study using 
Snodgrass' date, 10 Feb. 2009. The model, together with di- 
rect measurements of the tail within the camera's field of 
view, indicate it is constituted of dust grains with radii a 
ranging between a few tenths of millimetres to 20 mm, with 
a number of particles proportional to ^-3-44±o.08 (differential 
size distribution index). The total dust content of the tail is 
estimated to at least 8 x 10^ kg, which could be packed in 
a sphere with a radius of 40 m (assuming the same density 
p = 3 000 kg m"^ as above). 

- The position and shape of the tail's envelope was modeled as 
originating from dust emitted on a cone with a half-opening 
angle a ~ 40° pointing below the orbital plane, toward the 
forward direction, with a ejection velocity Ve < 0.20 m s"^. 
The X-shaped arcs were modeled as a collection of large par- 
ticles, with Ve > 0.20 m s"^; they were independently found 
to be emitted on a section of the same hollow cone. The are 
seen as an arc as a consequence of the geometric foreshort- 
ening of the cone section (after one year of evolution); other 
sections of the cone, which are not foreshortened in the same 
way, contribute to the overall tail. The Rosetta observations 
(ISnodgrass et al.ll20TQl) did not see other sections of the cone 
as arcs because of the coarser resolution of its camera. 

- Because ice sublimation, electrostatic lifting and rotational 
spin up were rejected as possible causes for the dust emis- 
sion, and because of the very short release time and the ge- 
ometry of the dust emission on a hollow cone with an half- 
opening angle a = 40°, we conclude that the dust release 
was caused by an impact by a small object. 

- Considering the volume of ejected dust in the framework 
of a gravity-dominated crater formation process, a crater of 
~ 100 m radius was produced in a total time of the order 
of seconds. While that crater is large compared to the body, 
other bodies show craters as large as themselves. An oblique 
impact — which is statistically the most likely — would ex- 
plain how P/2010 A2 escaped complete, catastrophic disrup- 
tion. 

In summary, the tail of P/2010 A2 can be explained by an 
impact on a ~ 80-90 m asteroid. The impact released dust in 
a hollow cone, which evolved for about one year to form the 
observed dust distribution. The fine-grained dust formed a long 
tail that was extended by solar radiation pressure, and the large 
grains (un-aff'ected by the solar radiation pressure) retained the 
conical shape. The viewing geometry revealed sections of that 
cone as the arcs in the head of the dust cloud. Therefore, what 
we witnessed was likely an event of collisional grinding within 
the Flora collisional family, and not a genuine, water-ice driven 
Main Belt Comet. With increasingly sensitive sky surveys cov- 
ering larger fractions of the sky during each dark period, we can 
expect to detect more and more of these collisions, as illustrated 
by the recent discovery of the dust cloud around (596) Scheila. 
While the original goal of this project was to probe the water ice 
content in the asteroid main belt, its end result was the observa- 
tion of a natural impact on a small asteroid. 
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